Resonant cancellation of off-resonant effects in a multilevel qubit 
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Off-resonant effects are a significant source of error in quantum computation. This paper presents 
a group theoretic proof that off-resonant transitions to the higher levels of a multilevel qubit can be 
completely prevented in principle. This result can be generalized to prevent unwanted transitions 
due to qubit-qubit interactions. A simple scheme exploiting dynamic pulse control techniques is 
presented that can cancel transitions to higher states to arbitrary accuracy. 
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Successful quantum computation depends on the ac- 
curate manipulation of the quantum states of the qubits 
|EJ. In practice, qubits are subject to many sources of 
quantum errors including thermal fluctuations of the en- 
vironment 0, qubit-qubit interactions ||, and intrinsic 
redundant degrees of freedom within a qubit such as the 
quasiparticle conduction in the superconducting qubits 
||,[| , and the effect of the higher levels in many practical 



qubit designs |3y£j] ■ This paper proposes a dynamic pulse 
control technique that efficiently eliminates unwanted off- 
resonance transitions. 

Various schemes to protect the qubit from qubit er- 
rors have been proposed that can be divided into two 
categories. The first one is the quantum error correc- 
tion technique ]t|-[T^] where the qubit state is encoded 
by redundant qubits. Different errors project the qubit- 
extra-qubit system into different subspaces that can be 
determined by measuring the state of the extra qubits. 
By applying a transformation according to the measure- 
ment, the correct qubit state can be restored. This ap- 
proach relies on large numbers of extra qubits to keep the 
errors from propagating. The second approach exploits 
'bang-bang' control techniques jl3| where the dynamics 
of the qubit and the environment is manipulated by fast 
pulses that flip the qubit. With the influence of the en- 
vironment averaged out, the qubit evolves in the error- 
free subspace. This method relies on the ability to apply 
pulses rapidly compared with the correlation time of the 
environment. This is an open loop control method. 

A particularly important form of intrinsic qubit error 
comes from the off-resonant transitions to the higher lev- 
els of a qubit when the qubit is being operated. Real 
qubits are not S = 1/2 spins that are perfect two level 
systems; redundant levels always exist that affect the 
information content of the qubit. The additional in- 
teraction that is introduced to achieve qubit operation 
by coupling the lowest two states of the qubit almost 
always includes unwanted couplings between the lowest 
two states and the higher levels. When the interaction 
is applied with frequency u> = L02 — u)% , resonant transi- 
tion occurs between the lowest two states; meanwhile off- 
resonant transitions to the higher states are also switched 



on. These transitions deviate the phase and amplitude 
of the qubit state from perfect Rabi oscillation. Numeri- 
cal simulations on the superconducting persistent current 
qubit (pc-qubit) |I|,[l4j show that this deviation can be se- 
vere when the unwanted couplings are of the same order 
as the Rabi frequency. 

In this letter, we study the effect of the higher levels 
on qubit dynamics during qubit operation by a group 
theoretic approach. We prove that the errors can be 
completely avoided by applying a time varying opera- 
tion Hamiltonian. Then we generalize this result to the 
qubit-qubit interaction problem which can be mapped ex- 
actly onto the first one. Using the idea of dynamic pulse 
control fl^l , we design a pulse sequence that cancels the 
leakage to the higher levels to arbitrary accuracy with 
0{N) number of pulses, N being the number of higher 
levels. 

The leakage to higher levels has two significant char- 
acteristics. First, unlike the environmental fluctuations 
that affect the qubit only slightly (less than 10~ 4 ) within 
one operation, the leakage changes the qubit dynamics on 
a time scale 1/wo that is much shorter than the qubit op- 
eration time (about I /w Rabi)- Conventional quantum er- 
ror correction corrects errors that occur with small prob- 
ability and is not a suitable strategy to cancel the fast off- 
resonant transitions. Neither can we use the bang-bang 
method to average out p| the leakage simply by manip- 
ulating the lowest two states, as the pulse induces these 
unwanted transitions at the same time. Second, ignor- 
ing all interactions with external variables, the leakage is 
coherent, although the coherent oscillation will collapse 
since the revival time is too long to be observed due to 
the large number of transitions of different frequencies 
]l6| ]. As will now be shown, the coherent nature of the 
leakage implies that this type of error can be corrected 
by applying a control sequence that coherently modifies 
the qubit dynamics. 

Consider a N level quantum system with the Hamil- 
tonian 7io, the lowest two states of which are chosen as 
the qubit states | f) and | |). The unitary transforma- 
tions on this N dimensional Hilbert space form the N 2 
dimensional compact Lie group U(iV). Without other in- 
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teraction, the trajectory of the qubit follows the Abelian 
subgroup {e~ lUot ,t G R}. 

Now apply to the qubit the perturbation "Hi, 
[Ho, Hi] =/= 0, to induce a desired transformation of 
the qubit. In most physical systems, unwanted tran- 
sitions to the higher levels are simultaneously induced. 
For example, in the pc-qubit Q operation, (Hj) mn = 
2nSf(m\ sin {2<p m + 27r/)|n) cos tot, when the bias flux is 
modulated with rf components of amplitude Sf and fre- 
quency uj. This perturbation has couplings between all 
the levels. By successive commutation of Ho, Hi, and 
their commutators until no independent operator ap- 
pears, a Lie algebra Ai is created. In almost all cases, 
Ai = u(N) @, u(N) being the Lie algebra of U(iV). 
The only exception occurs in a zero measure subspace of 
u(iV) when Hi and Ho are both in the same subalgebra of 
u(iV). Thus, with almost all perturbations, the evolution 
operator can be any element in \J(N); and transitions to 
higher levels are unavoidable with an initial state that 
only occupies the lowest two levels. 

To prevent the transitions to the higher states at time 
t means to restrict the evolution operator U(t) to the 
submanifold of U(2) © U(A - 2), U(2) being the uni- 
tary group on {| t), | J,)} and XJ(N — 2) on the remaining 
N — 2 states. This applies 4(N — 2) real domain re- 
strictions onU(t): U{t) lk ,U(t) 2k = 0,k = 3 ... N. In 
contrast to a perfect qubit operation during which U(t) 
remains in the subspace U(2) © XJ(N — 2) all the time, 
the qubit is allowed to stray away from this subspace if 
only it goes back to this subspace at the designated time 
t. The qubit dynamics can be manipulated by varying 
the strength and phase of the perturbation with time. 
As the N 2 dimensional Lie group U(iV) is compact, any 
transformation can be reached at time t by adjusting the 
N 2 + 1 parameters in the following process p7| : 

where a is introduced to ensure that t = Y)tj. By play- 
ing with the N 2 + 1 real parameters, the 4(7V — 2) real 
numbers in U{t)ik,U(t)2k can be set to zero so that the 
state of the qubit stays in the {t, 1} space without leak- 
age. Hence by turning the perturbation on and off 0(N 2 ) 
times, the lowest two states are completely decoupled 
from the higher states. 0(N 2 ) pulses give a sufficient 
condition that is required to achieve arbitrary transfor- 
mation. As will be shown later in this letter, with proper 
arrangement, we can design a pulse sequence of O(N) 
pulses to cancel the transitions to the higher levels. Un- 
like that in the quantum Zeno effect Jl8| where measure- 
ment is used to prevent the system from evolving, the 
dynamics in this process is described completely by uni- 
tary evolutions. 

As the unwanted transitions are off-resonant transi- 
tions whose amplitudes decrease roughly as 7^ / u>i (uii is 
the energy of the iih level, 7^ the coupling between level 
i and j), the influence of the levels with LOi/bJo 1 can 
be ignored. In the pc-qubit M the energies of the lower 



levels increase fast enough (ljxo > lOwo) that levels be- 
yond 1 10) can be ignored. The energy of the ith level of 
the charge state qubit || increases as i 2 ; less levels affect 
the qubit dynamics than that in the pc-qubit. Hence the 
number N of the higher states that are involved in the 
qubit dynamics in real designs can be reasonably small. 
As a result, the number of pulses in the previous analysis 
is also reasonable. 

One question to ask is whether there is any fundamen- 
tal difference between the errors due to transitions to the 
higher levels and those due to the fluctuations of the en- 
vironmental variables. Putting it in another way: what 
is the difference between the intra-qubit coupling in a 
multilevel qubit and the qubit-external-system coupling? 
In the following we will show that the iV-level qubit can 
be mapped into interacting subsystems, and vice versa. 

Let the initial state of a N- level qubit be |^o) = 

|1) + c4°' ) |2), |1) an( i 1 2) being the lowest two states. 
To map the qubit into two subsystems, we divide the N 
states into two subspaces SP\ and SP 2 by adding the 
vacuum states \V\) and | V2) to the respective subspaces 
as SPi = {\Vi), |1>, |2)} and SP 2 = {\V 2 ), |3), . . . , \N)}. 
Now the N dimensional Hilbert space of the original 
qubit is embedded in the 3(N — 1) dimensional direct 
product space SP± © SP 2 . The states in the expanded 

space are \W) = fcM^fif ')» where ^ and h f 
are basis of the two subspaces respectively. The initial 

state is |* > = (ai 0) |l) + a ( 2 ] \2))\V 2 ) in the expanded 
form. The unitary transformations in this expanded 
space forms the group U(3(iV — 1)). 

Perturbation introduces coupling between different 
states. When mapped to the expanded space, the per- 
turbation Hi connects states in the N dimensional sub- 
space spanned by {\1)\V 2 ), \2)\V 2 ), |Ui)|3), ... , \Vi)\N)}. 
So Hi and Ho create N 2 dimensional subalgebra u(iV) 
in the expanded space. Under the perturbation, the 
wave function in the expanded space can be described 
as \W) = ( ai |l) + a 2 \2))\V 2 ) + Y%=3 <*i\Vi)\i) , where a, 
are time dependent parameters evolving with the pertur- 
bation. 

From this analysis, the higher levels in the qubit form 
an effective environment that interferes strongly with 
the lowest two levels. Interaction strength is the major 
difference between this effective environment and a real 
one The couplings between SP\ and SP 2 are strong 
and comparable to the Rabi coupling that realizes qubit 
operation. In contrast, the interactions between the en- 
vironmental oscillators and the qubit are weak due to the 
0(l/yV) factor that originates from the normalization 
of the extented modes || . So the thermal fluctuations are 
not enslaved to the qubit dynamics and can be treated 
classically. The strong interaction with the higher lev- 
els also explains why the error due to leakage occurs at 
such a short time that a particular strategy is required to 
correct the error. Another thing to mention is that this 
effective environment only comes with qubit operation, 
while the real environment affects the qubit all the time. 
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Hence we worry about the leakage only during qubit op- 
eration and choose to correct the leakage by controlling 
the operation process. 

By reversing the mapping, interacting qubits can be 
modeled as one multilevel quantum system. One example 
is two interacting qubits with basis jb^ ),« = 1 . .. Ni, 

(2) 

and | bj },j — 1 ... N2, respectively. Labeling the state 

\bP)\bf } ) as \(i - l)N 2 + j), we have the Aq7V 2 basis 
for the equivalent multilevel qubit. In the same way, 
n two-level qubits form a quantum system of 2™ levels. 
Here the number of states grows exponentially with the 
number of qubits as the entanglement between qubits 
has to be included in a single system now pjj . Pertur- 
bation applied to one of the qubits can cause unwanted 
couplings within the 2™ levels, and induce off-resonant 
transitions that affect the qubit performance. Similar to 
the couplings of the multilevel qubit, these couplings are 
also strong and cause fast errors. Taking the pc-qubit 
as an example, the interaction between the two qubits is 

§: H 2 = J,<t£ ) (t { ? ) + Jxial^a™ + a^ajP), where J z 
and J x terms are due to the inductance coupling between 

qubit circuits. When a term is applied to rotate the 
first qubit, the second qubit will be involved and qubit 
dynamics will be changed. 

Although the mapping of the multilevel qubit and the 
multiqubit system into each other is just another way 
of looking at the same problem, it shows that the er- 
rors due to the qubit-qubit interactions || can be treated 
by the same approach as is used in cancelling the inter- 
ference of the higher levels. Again we turn to the idea 
of coherent pulse control that is exploited in the higher 
state problem. Now the number of pulses increases ex- 
ponentially with the number of interacting qubits, but it 
doesn't cause a disaster in real designs where only the 
nearest neighbour qubit interactions are important and 
n can be made small in the qubit layout geometry. 

To illustrate the general idea of exploiting dynamic 
pulse control to cancel the errors due to the higher states, 
we give an example of pulse sequence that completely 
cancels the transitions to the higher levels with 0(N) 
pulses. Let us start from a three level system with eigen- 
values u>i, i — 1,2, 3. The energy difference between level 
i and j is shorthanded as LOij. An interaction 7ij that 
couples level i and j by 7^ is applied to operate the 
qubit. When the third level is not present, 712 is the 
Rabi frequency of the lowest two states. For simplicity, 
we ignore the diagonal couplings 7^ as Cw t . As will 
become clear, the effectiveness of the designed pulse se- 
quence depends on the condition \jij/u>ij\ <SC 1 which is 
satisfied in most qubits. 

The Hamiltonian in the interaction picture is Hint = 
e lHot Ttie~ lHot cos (uit + (/)), uj being the pulse frequency. 
The wave function ^(t) = [uvw] T evolves according to 
the equation = Hi n t^{t)- When the perturbation 

is weak, this equation is integrated order by order as: 



= *(0)+ f*dt"H int (t')y(0) 

+ /* dt' dt"H m t(t')H mt (t")*(0) 



(2) 



The cosine function is used in the rf pulse instead of 
the single frequency rotating wave. In many systems, no 
physical correspondence of the rotating frame exists. For 
example, the circuit of the pc-qubit is biased by z direc- 
tion magnetic flux and the perturbation is high frequency 
modulation of the z flux. No rotating frame of transverse 
field can be defined for the oscillating flux. 

Our strategy to reduce the unwanted transitions is to 
divide the qubit operation into short intervals of t and 
attach additional pulses to each operation pulse to cor- 
rect errors from this short interval. The operation pulse 
is in resonant with CJ21 of the lowest two states. Be- 
sides rotating the qubit between the level 1 and 2, it 
brings up off-resonant transitions between the third level 
and these two levels through the couplings 713 and 723. 
Then the same perturbation is applied in two other pulses 
with different frequencies, amplitudes and phases as: 
a^iHi cos (w 3 it + 03i ) and a&Hi cos (w 32 t + ^32), both 
for time to cancel the unwanted transitions to the 
third level. This three-piece sequence is repeated T op /to 
times to finish the qubit operation. The time to satisfies 
1/lo <C to *C — 1,2,3 with both l/w2ito and 

"fijto being small parameters of the same order. Thus we 
have two small parameters in this process. This is crucial 
for this simple pulse sequence to work. 

Starting with an initial wave function 5'(0) = 
[uqVoWo] t , wo — 0, after the CJ21 pulse, the third level 
has the component: 

W = Unf" ^ ; '-) 

7 » 3(e -i^2i--32)t0-i) ^(e'^l+^O,!) ( (3) 
1 1 ■ OJo-i — OJio 1^21+^32 



u Q 9 u + v 9 v 



where 6 U and 9 V are of third order. The main compo- 
nents in w are second order terms that depend on the 
initial wave function uo and vq linearly. With to satisfy- 
ing e 2lW21 *° = 1, u and v have third order deviations from 
the correct two-level rotation. The other two pulses are 
then applied to cancel the w component. The w 3 i pulse 
induces a resonant transition between level one and level 
three to cancel the uq term in w; the w 3 2 pulse induces 
a resonant transition between level two and level three 
to cancel the ^o term in w. The amplitudes and phase 
shifts of these two pulses can be expanded in ascending 
order as: 



'32 



*31 

(2) i<A (2) 
*32 e 



(4) 



The first order coefficients cancel the second order terms 
in w and modify the higher order terms 9 U and 9 V when: 



(1) i<j>^ _ e -'("21""3l)«0_ 1 _ e i("21+"3l)'0 



(1) i<A (1) 

% e 32 



i(a)21— "3l)*0 
i(u> 2 l -^32)*0 _1 

i(ld21— ^32)to 



i(w21+W3l)*0 
e i(w 2 l+^32)*0_l 

i(u21+<*>32)*0 



(5) 
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It turns out that the nth order terms of w after the 
two correction pulses include linear terms of cv 3 ™ and 



«32 ^ , and complicated terms that depend on a 



MM 



3i 

(k = 1 . . . (n — 2)). So, for any n, aJj™ 1 and a 3 2 
can be determined by the lower order components of 0:31 
and 0:32 to cancel the nth order of w. As a result, the 
transitions to the third level can be completely erased. 
The parameters a^i and 0:32 do not depend on the ini- 
tial wave function uq and vq- This is similar to solving 
the wave function in the perturbation theory where the 
higher order terms are derived after the lower order ones. 

After the fcth pulse sequences, with w = 0, the wave 
function of u and v is: 



Uk+l 




cos ip + s u —ism(p + t u 




Uk 


. Vk + 1 _ 




— i sirup + i v cos <p + s v 




Vk 



(6) 



where tp = 712^0 is the phase rotation of the two- 
level qubit; the s and i terms are of third order. As 
w = 0, this is a unitary transformation that devi- 
ates from the original Rabi oscillation by third order 
corrections. The matrix can be written as l4(to) — 
exp (— i(ji2cr x + Sq + JV fl»Oi)io), where Si are third or- 
der small numbers that can be determined by known pa- 
rameters and do not depend on the index k. This is a 
renormalization of the qubit operation 7 12 with the third 
level decoupled. 

This correction strategy is easily generalized to N(N > 
3) level system. By applying rf pulses with frequencies 
0Jii,uii2,i — 3...N, the transitions to the higher levels 
are completely erased. Assuming no particular symme- 
try between the states, 2 (TV — 2) pulses are required in 
this process. 

One may wonder why this simple pulse sequence works 
so well to correct the transitions to the higher states. 
For N — 2 higher levels, to decouple these levels is to 
exert 4 (AT — 2) real domain restrictions on the transfor- 
mation matrix: Un,Un = 0, i = 3...N. Our tools are 
the Hamiltonians Ho and Tii that create the whole u(N) 
algebra by commutation. Our pulse sequence U(to) = 
n 4i/3 P(a l/3 ,^ /3 )e-^^ cos ^ 1 *' dt ' (i = 3 ... N and =| 
,| ), P(a t p,<j> ll3 ) = e -i/«/«i/3co S (c^t'+& /3 )dt' ) contents 
A(N — 2) free parameters. By choosing proper pulse se- 
quence, we can achieve the decoupling with proper pulse 
parameters. 

In conclusion, we discussed the errors due to unwanted 
transitions to the higher states of a qubit during qubit 
operation. It was shown by a group theoretic argument 
that these errors can be completely prevented in princi- 
ple. Then we generalized the result to the errors due to 
qubit interactions, which can also be prevented when the 
number of coupled qubits is not large. A simple pulse 
sequence that modifies the qubit dynamics and cancels 
off- resonant transitions to arbitrary accuracy with 0{N) 
pulses was proposed to illustrate the general analysis. 
Our results showed that the idea of dynamic pulse con- 
trol MM also works for the fast quantum errors due to the 



higher states of a qubit. These results suggest that dy- 
namic pulse control, together with conventional quantum 
error correction, can function as a powerful tool for per- 
forming accurate quantum computation in the presence 
of errors. 
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